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Abstract 

The well-established effective action and effective potential frame- 
work from the quantum field theory domain is adapted and success- 
fully applied to classical field theories of the Doi and Peliti type 
for diffusion controlled reactions. Through a number of benchmark 
examples, we show that the direct calculation of the effective po- 
tential in fixed space dimension d = 2 to one-loop order reduces to 
a small set of simple elementary functions, irrespective of the mi- 
croscopic details of the specific model. Thus the technique, which 
allows one to obtain with little additional effort, the potentials for 
a wide variety of different models, represents an important alterna- 
tive to the standard model dependent diagram-based calculations. 
The renormalized effective potential, effective equations of motion 
and the associated renormalization group equations are computed in 
d — 2 spatial dimensions for a number of single species field theories 
of increasing complexity. 
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1 Introduction 

The effective potential provides physical insight into the behavior of non- 
equilibrium stochastic dynamics, and takes into account both nonlinear 
interactions and random fluctuations to a given number of loops Pp. It 
has been used to study transitions in non-equilibrium steady states |2 by 
means of the renormalization group improved equation of state T , and has 
been used most recently in nonperturbative renormalization group (RG) 
methods applied to nonequilibrium critical phenomena 0]. 

Given the importance of the effective potential in statistical physics, it 
is useful to have efficient analytic calculational methods for its evaluation 
that are applicable to a wide variety of non-equilibrium field theories. In 
this paper, we take the well-established effective action and effective po- 
tential formalism from quantum field theory and adapt it to obtain renor- 
malized closed form analytic expressions for the perturbative potential and 
attendant renormalization group equations (RGE). Once this is carried in 
detail for one given model, it is easy to extend it to a wide class of different 
reaction models with little additional effort. 

To illustrate and validate these techniques, we will consider some well 
studied field theories of reaction and diffusion processes that have been 
derived from the microscopic master equations following the procedures of 
Doi [H] and Peliti 6 . We will show by explicit example that the method 
allows one to obtain both renormalized potentials and RGEs with little 
effort. In fixed dimension, this represents a major advantage over the deli- 
cate diagram based approach where each separate reaction model requires 
defining and handling a distinct set of diagrams and the working out of 
combinatorics. 

We now outline the basic steps that take us from a given bare action to 
the one-loop effective potential. The field theories we are interested in are 
characterized by bare (unrenormalized) actions S[(j), (j)] that depend on a set 
of fields (jji, an equal number of response fields where i = (1,2, ...,iV) 
and N is the number of independent species in the model. From S we 
construct the following array of second order functional derivatives of the 
action: 



where i,j — (1, 2, ... , N). This is written in 2 x 2 block form. Since the 
potential is defined for constant and steady field configurations, we easily 




(1) 
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evaluate the Fourier transform (FT) of Eq. ^ via 

\ <f).<f)— const./ 

(2) 

where M is a 2N x 2N matrix and d is the spatial dimension. We then 
evaluate the determinant of M and insert it into the expression for the 
one-loop potential (see the Appendix for a self-contained derivation): 

f r 1 ^ detM(fc,u;;0,0) ^ 

V(0, 0) -Vo{c^,c^) + jJ y_ ^ In ( detM(fc,c.;0,l) J ' 

Thus, given an action S[4>,(j)\, it is straightforward to calculate the cor- 
responding matrix M, Eq. |(2Jl, followed by its determinant, and thence 
to the effective potential Eq. ((2Jl. The zero loop contribution Vo corre- 
sponds to the non-derivative part of the action S (see Eq. (|133ll '). When 
this is renormalized, we will refer to it as the tree-level potential. The one 
loop contribution is proportional to A. The indicated integrations over 
frequency u and wavevector k can be often be expressed entirely in terms 
of elementary functions for fixed dimension. These points will be demon- 
strated explicitly in what follows. 

Once the effective potential is calculated to some order in perturbation 
theory or loops, it is a simple matter to obtain the associated RGEs [7]- 
This procedure for obtaining the RGE's is based on the so-called derivative 
technique, which was first introduced in the work of Fujimoto et al. [5] and 
subsequently extensively developed in The important observation is 
that the bare theory does not depend on the arbitrary scale /x introduced 
by the renormalization scheme. This observation leads to a simple identity. 
This identity yields an equation which is a polynomial in the background 
field whose coefficients are precisely the RGE's that we seek. This is a very 
simple and powerful technique that deserves wider use. Equally notewor- 
thy is the fact that to one-loop order, the case we treat here, the analytic 
calculation of the effective potential is entirely straightforward and one can 
dispense with model-dependent diagrammatic manipulations. This has the 
clear advantage of allowing one to survey and analyze comparatively a wide 
variety of reaction-diffusion models by making simple algebraic substitu- 
tions to some of the terms appearing within the universal formula for the 
one-loop effective potential. Knowledge of the effective potential also allows 
one to reconstruct the effective action (up to wavefunction renormalization) 
and deduce the corresponding equations of motion obeyed by the physical 
density field. These equations of motion are deterministic, yet incorporate 
the presence of the internal reaction noise through modified effective force 
terms. 

The remainder of this paper is concerned with the explicit evaluation 
of Eq. (PJ for a selected number of single species [N = 1) diffusion limited 
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reactions in d = 2, and is organized as follows. In the next Section, we 
calculate the one-loop effective potential corresponding to pure pair anni- 
hilation A + A^(lj. This is renormalized and the resultant finite expression 
is examined with and without the loop corrections to reveal qualitatively 
the role of fluctuations on the asymptotic states of the system. With the 
potential calculated, it is a simple matter to obtain the associated renor- 
malization group equation (RGE) that controls the scale dependence of the 
annihilation reaction rate. We also derive the effective equation of motion 
satisfied by the fluctuation averaged particle density. In Sec. 01 we calcu- 
late the effective potential for the Gribov process which is described by the 
set of single species reactions A^(ll,A^A + A, A + A^A and extract 
the one-loop RGB's and derive the effective equation of motion. Reaction 
and diffusion systems with the underlying reaction process A + A ^ (d and 
A —^ (m -|- 1)A are known as branching and annihilating random walks 
(BARW). Here we calculate the one-loop effective potential for arbitrary 
even m in Sec. ^and extract the RGB's in Sec. 14.21 We then specialize to 
the case m = 2, solve the corresponding RGB's and obtain the equation of 
motion for the noise-averaged density field. The purpose of this paper is to 
carry these steps out explicitly in detail in order to illustrate the techniques 
involved and the ease with which they are applied to a variety of reaction 
models. Conclusions and discussion are presented in Sec. |S1 

2 Pair annihilation reaction A + A ^ ^ 

Perhaps the simplest of all two-body irreversible recombination processes is 
the single species pair annihilation reaction A + A inert. In this model, 
density fluctuations were shown to alter radically the asymptotic decay laws 
predicted from the mean-field rate equation approach for all dimensions 
d < 2 through a series of dimensional considerations, heuristic and scaling 
arguments and computer simulations ^01^3- The anomalous scaling in 
the density was firmly established by Peliti on the basis of field theoretical 
descriptions using the renormalization group fRG)|12|. Subsequent works 
applying the RG to this system include El ■ An RG calculation 

for the reaction kA inert is carried out in jJB], which recovers Peliti's 
result for k ~ 2. To illustrate the techniques mentioned above, we present a 
detailed calculation of the one-loop approximation to the effective potential 
in this simplest possible case. We renormalize and examine this potential, 
use it for a simple derivation of the corresponding RGB's in d = 2, and 
to obtain the effective equation of motion for the noise-averaged particle 
density. 
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2.1 One loop effective potential 

The following bare action will be the starting point for our effective poten- 
tial calculation: 

S[tp,4>]:^ Jd'^xdt ^(dt+^-D\/^^tl>-X{l-i:^)i;^ . (4) 

The continuum parameters D, A appearing in the action Eq. (0J stand for 
particle diffusion and the annihilation reaction rate, respectively. A field- 
theoretic action S corresponding to the pair annihilation reaction including 
diffusion was first written down almost twenty years ago by Peliti |12| , but 
without the inclusion of the bare parameter ^. This deserves some special 
mention. Physically, this is a mass term, and is normally taken identically 
to zero from the outset in field-theoretic treatments of pair-annihilation. 
Here, we must explicitly add it to the bare tree-level action Eq. since, as 
we will demonstrate below, it is required in order to carry out consistently 
the one-loop renormalization of the effective potential. Nevertheless, after 
this renormalization is performed, we are free to set the finite and renor- 
malized value of ^ to zero, and we will do so. Our renormalized action 
agrees with the renormalized action obtained by other means ^2 • 

We next carry out the basic steps outlined above in Sec (Q) that take 
us from the bare action to the formal expression of the one-loop potential. 
The matrix Eq. of second functional derivatives of the action is 

2X(f>'^ _ ~dt+^- DW'^_+ 4A0(/) 
dt+^- DV^ + AX4>(t> -2A(1 - 

X 5'^{x-y)&{t- s). (5) 



This is a tensor product of a 2 x 2 matrix in field space times an infinite 
dimensional but diagonal array in configuration space and in time (the 
product of delta functions). Since the effective action and potential depend 
on the averages of the fiuctuating fields, we must distinguish between the 
pairs of fields <f), ip, and (l),-ip, respectively; see Eq. (I128f) . For constant fields 
we pass to Fourier variables to obtain the array M, Eq. Q, where 

2A02 _ ^iuj + ^ + Dk'^_+ 4A00 

^ + DP + AX4>(j) -2A(l-(^2) 

(6) 

Inserting this matrix into the expression for the one- loop effective potential 
Eq. ||2Jl yields (recall that Vq — yt ^'-"^ constant fields) 

2/1 22 



M(fc,^;0,^)= ^,,,^,^^,2 , o^n 12, 



Vi^,cj))^-X<j,\l-^') + 

A f d<^k /"^ + [i:)fc2+^ + 4^202(1 _ 02) 



(7) 
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We first carry out the frequency integral, wiiicii can be evaluated exactly 
in closed form. By means of the following identity (see Eq. (4.222.1) of \17\) 
valid for C,A> 



cLu In 



we can write the one-loop 0{A) contribution in Eq. (TJ as follows: 



^ " 2X fc'^-^dfcy ^ |^(m2+^ + 4A#)2 + 4AV(l-02) 
- (9) 

where dfld is the element of solid angle in d-dimensional wavevector space. 
We now temporarily regulate the integral over wavenumber modulus k with 
an ultraviolet cut-off fc < A and evaluate it directly in d = 2 dimensions. 
Introducing the change of variable u — Dk^, we can express the regulated 
d = 2 integral Eq. as follows: 



^ ru=DA^ 



8nD 



du 



+ am + aa - \/iu + O^j , (10) 



where 



ai = 2(^ + 4A#) (11) 

a2 = ^-f4AV(l-0^) (12) 

A = 4a2 - a? = 16A2,/,2(1 - 02^, (13) 

and which can be expressed in terms of the following elementary functions 
(see e.g., (2.262.1) and (2.261) in [n|): 

{87tD)I^ = «i) ^ ^^^^2 ^ (j^A2)2 (14) 
+ — In (2y^a2 + aiDK"^ + {DK'^f + 2DK^ + ai 



- i(i^A2 + e)^-|v^+^-|ln(2V^+aO. 

A careful asymptotic analysis of the regulated loop integral Eq. (|14|l reveals 
that Ia contains a quadratic and a logarithmic divergence in A in addition 
to certain finite terms as this UV cutoff is taken to infinity, namely 



- ^V^}+0(A-2). (15) 
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A particularly useful separation between finite and divergent pieces is ob- 
tained by introducing an arbitrary finite momentum scale fj, and writing 
the argument of the logarithm as follows: 

After using this simple identity Eq. (|16|l . we have the following regulated 
one loop effective potential Eq. Q in d = 2 space dimensions: 

v(</>,0) = -Xc^\i-4>') + 

A 



SttD 

+ -L(A + 2a?)-^V^|+0(A-2). (17) 

This expression is given in terms of bare (unrenormalized) parameters (A, ^) 
and the dependence on the UV cutoff A is explicit. As pointed out above, 
the tree-level term S,4'4' will be needed for the renormalization program. 
Note that the regulated effective potential Eq. ((T7jl . is seen by inspection 
(due to elementary properties of the logarithm function) to be manifestly 
independent of the arbitrary scale fj,. This is the key observation needed in 
order to derive the RGB's, of which we will make direct use of below. 

We now proceed to renormalize, i.e., absorb the divergences into the 
bare parameters A and ^. To carry this out, we write 

A = + Aki + 0{A^) (18) 

e = a^^) + AK2 + 0{A^), (19) 

in which we decompose each bare parameter into a finite renormalized and 
possibly /L(-dependent part plus a divergent counterterm. Note the divergent 
counterterms start off at order 0{A), which is the same order at which the 
divergences in the regulated potential appear. These counterterms ki,K2 
are chosen to render the one-loop effective potential finite and UV-cutoff 
independent. Note the very important fact that the stochastic fluctuations 
induce divergent terms at the one-loop level whose coefficients have iden- 
tical polynomial field structure (namely, ^ — <^ and A, see Eq. (|ll|l and 
Eq. H13(l l as those appearing in the bare tree-level potential Vq. Inserting 
Eq. I|18I1H(I into Eq. ^ yields the solut ion for the two counterterms: 

« = -4|^A^ (21) 
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In this final step, and after cancelling off the UV divergences, we obtain the 
finite and renormalized one-loop effective potential for pair-annihilation (in 
d= 2): 



= -A(/i)0'(l-0') 




(22) 



We are now free to assign zero to the finite renormalized value of ^, and 
have done so here and in the following subsections. Below, when we come to 
discuss the RGB's, we will see that the renormalized parameter ^ does not 
run with scale, at least to one-loop order, proving that this is a valid choice. 
In this expression Eq. H22|) it is understood that A = A(/i) is the renormal- 
ized, /Lt-dependent coupling and we have set the finite renormalized value 
of C(/x) identically to zero. We next determine their scale dependences. 



Since the bare theory Eq. (|17|l does not depend on the arbitrary scale 
/X introduced by the renormalization scheme (recall our use of the trivial 
identity in Eq. (|15|l ). we are able to derive the renormalization group 
equation that the coupling A(/i) must satisfy, from the identity [7| HJ |H] 



which upon inserting Eq. (|22() immediately yields the RGE. Differentiating 
as indicated, we get 



2.2 The one-loop RGEs 



rfV(0>) 



= 



(23) 



^ .|-(-A^^(l-^^)) 




4- X^<j>^{l - (f?) + 8A2(^^^2 _ 4A2(/,2^y^ I + 302^ + 0{A^), (24) 



which implies 



A 



\'(l)''{l-4>'') + 0{A^). 



(25) 



2itD 
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Note that derivatives with respect to ^ of the other terms on the right hand 
side of Eq. proportional to A lead to contributions at the next higher 
order 0{A^), and so can be neglected at the one-loop order we are working 
at. Indeed, from Eq. H18|) we see that dX{^)/dii — d{X — Ani{^)) / dji = 
-AdKi{pL)/dn = 0{A^). 

Cancelling the common field-dependent factor — (j)^) from both 
sides of Eq. ^ yields the RGE for A in d = 2: 

^d^I ^ 2^^ ' 

valid to one-loop order. This can be solved immediately and gives 



A(m) = Ao(l-^lnii + 0(^2 
\ Zt:U hq 



where Aq = X{po), is the value of the coupling at some initial arbitrary 
reference scale fig. 

If we retain the dependence on the renormalized ^ in Eq. (|22|l and de- 
duce the consequences following from Eq. (|23|l , we now obtain two RGE's, 
corresponding to the two independent polynomial terms (/^{l — c/r) and (j)(j>, 
respectively. Once again, we obtain Eq. (|26l) as well as 

= + O{A^), (28) 

c(m) - m + oiA^), 

indicating that the parameter ^ does not run with scale, at least to one- 
loop order. So, it is consistent to take the renormalized (but not the bare) 
mass parameter to be strictly constant, and setting ^(/i) = ^(0) = is a 
perfectly valid choice. 

This RGE Eq. 1)26(1 has been calculated directly in d = 2, which coin- 
cides of course with the upper critical dimension of the pair-annihilation 
process ^Jj. When working with the RG, it is useful to have knowledge 
of the scale dependence in variable dimension d. This would require us to 
go back and attempt to calculate the potential also in arbitrary d, which 
is a difficult task. This is because the one loop potential represents an 
infinite sum of Feynman diagrams. While individual diagrams are usually 
easy to analytically continue to arbitrary dimensions, their infinite sums 
are not. However, we can gain some useful insight into the structure that 
the variable-d RGE must have by using some simple dimensional analysis 
in conjuntion with the d = 2 RGE in Eq. ]2<j\i . In this way, we arrive 
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at the approximate mathematical form of the RGE for the dimensionless 
coupHng, which may be expected to hold for values of d close to d = 2. 
We comment below in what sense it is approximate. In regards the scaling 
dimensions of the field theory, the diffusion constant D can be absorbed 
into a rescaling of the time t Dt. Doing so in the action Eq. Q and 
expressing the dimensions of the resultant terms in the action in powers of 
momentum yields jl6 | 

[t]^p-' m=p' m^p" [x/D]^p'-''. 

(29) 

From the last equality in Eq. (|29|l . we are lead to define a dimensionless 
coupling g in any d via 

^'W' (=•») 

where e — 2 — d. 

We can write Eq. H26|l so that it is dimensionally consistent for all 
dimensions d. This is achieved by first noting that (we set ^ = 1 in what 
follows) 

^d\ ^ / Kd 

djl \Dp'^)d=2 

where 

= J-^^-^^w^v^'^'-r' (32) 

is a geometric factor that results from the integration over solid angles and 
which is part of the complete d-dimensional loop integral in Eq. Q . Note 
that K2 = ^ and that g{d = 2) = ^. We next extend Eq. to all 

dimensions by writing 

n— = — — A . 33) 

^dn Dn'' ^ ' 

This equation is of course only true in d = 2, and the approximation 
comes in using dimensional analysis to extend it to any dimension. Lastly, 
substituting Eq. l|5n|l into Eq. (|^ . and expressing A in terms of 5, yields 
the (admittedly approximate) one loop beta function for g: 

/3(5) = Ai^ = -e.g + i^d.g2. (34) 

This beta function Eq. (|34|l is quadratic in g and for d < 2 has the nontrivial 
fixed point /3(5'*) = at 5* = = 0(e). This fixed point is of order e. 
For d > 2, = 0. Inspection of Eq. shows that the calculation of 
the one-loop RGE for A valid for arbitrary d would require the integration 
of the corresponding integrand over the modulus of k with weight k'^~^, 
in addition to the integration over solid angle Eq. H32(l . The above result 
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takes the latter in to account, but the former is evaluated at fixed d = 2, 
so we can expect the overall numerical coefficient of the quadratic term 
to be in error in this respect. Nevertheless, for values of d suflficiently 
close to d = 2, this should provide qualitatively correct results. By way of 
comparison, the exact one-loop beta function calculated in is given by 

,35, 

For d > 2, = 0. For d < 2 the nontrivial fixed point is given by 
9l = T{2-d/2) ■ At d 2, these two beta functions Eq. and Eq. (jSU 
are identical, and we have I3{g) — ^g^- For d < 2 we can compare the 

deviation of the approximate to exact nontrivial fixed point via the ratio 

* / * _ r(2-rf/2) 

2.3 Effective equation of motion 

From the one-loop potential we can construct the effective action to the 
same (one-loop) order via [TH| [El EH EH 122 



n 



= J d'^x dt 1^1/2^1/20(9* - L'V^)^ + V(<?!>(a;, i), <?!)(x, t)) -f . . .} , 

(36) 

where Z and Z denote the wavefunction renormalization constants for 4> 
and 0, respectively, and the dots represent higher derivative terms in the 
fields. From the fact that there is no wavefunction renormalization for 
pair- annihilation B*, '16', we can immediately set Z — Z — 1 in Eq. 
Then, from Eq. (|13()|l and Eq. H22|l . we can derive the (one-loop) effective 
equations of motion for pair-annihilation in d = 2 obeyed by the noise- 
averaged fields 4>{x,t) = (?/;(a;, t))(Q Q) and (f>{x,t) — {ip{x,t)) (q q^ We find 
that the equation of motion for (j) (recall = 1) 

is identically solved by the stationary and homogeneous solution (f> — 1. 
The equation of motion of the physical density field (p follows from 



V 56 



^ = 0, (38) 



which when written out explicitly, implies that 

i9t-DV^)^ = -S. , (39) 
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= F^. 

Here, the effective force experienced by the field is denoted hy F^ — —~. 
In Fig. ^ we show a quahtative example of the force experienced by the 
density field cj). Notice that since the recovery force is less strong in the 
renormalized case (the fluctuations cause the slope at the origin to be 
less negative, relative to the tree- level case), the rate of approach to the 
origin will be slower. Putting ^ = on the right hand side of Eq. H39(l 
yields the so-called rate equation which one might write down as a first 
approximation to the equation of motion for the density. If we let n{t) be 
the spatial average of (j), then we see that the only stationary solution of 
Eq. (|39|) is n = 0, corresponding to the inactive state with no particles 
remaining, which of course, comes as no surprise. Zeroes of the effective 
force correspond to maxima/minima of the effective potential. Therefore, 
since the logarithm of a small number is negative, it appears as though 
the one-loop correction in Eq. (|39|l has turned the minimum at the origin 
into a maximum and caused a new minimum to appear away from the 
origin, corresponding to a putative noise induced active state. However, it 
is wise to proceed with caution in this specific example: the apparent new 
minimum occurs at a value of n determined by (note from Eq. (|29|l that 
the ratios X/D and are pure dimensionless numbers in d — 2) 

Ax ,4An 
dJ \Dfi 



Since we expect higher orders in perturbation theory to bring in higher 
powers of {X/D) In {AXn/D^j?) , this new state hes well outside the expected 
range of validity of the one-loop approximation I22j . even for arbitrarily 
small reaction rates (or large diffusions) X/D. This new solution Eq. H4()|l 
must be rejected as an artefact of the approximation ^. Higher loops are 
needed to assess the situation. However, this same criticism can not be 
applied to the effective force associated with the Gribov process, the next 
example to be considered below. 



3 Gribov process ^ ^ 0,^ ^ ^ + ^,^ + ^ ^ ^ 

As a second application of the methods outlined in Sec we consider 
the competing reactions of single particle annihilation, splitting and two- 

^Tho one- loop effective force in the pair-annihilation model is mathematically similar 
to the one-loop potential for the massless quartically self-interacting meson field, for 
which an analogous criticism applies. This point has been taken over and adapted to 
the context of the present paper; see Ref. I22| . 
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body recombination. This is a simple model exhibiting a non-equilibrium 
phase transition between active and inactive or absorbing stationary states 
in which the stochastic fluctuations cease entirely. A mathematical equiv- 
alence between such chemical processes and Reggeon field theory has been 
firmly established |^ |^ . Below we calculate the corresponding one- loop 
potential, carry out the renormalization and obtain and solve the RGE's 
for the rates of pair recombination and fission directly in d = 2 dimensions. 

3.1 One loop effective potential 

The following field-theoretic bare (unrcnormalizcd) action S will be the 
starting point for our effective potential calculation: 

- X{l - 11!)^;^^ + biP + a] . (41) 

The continuum parameters D, X appearing in this action stand for particle 
diffusion and the reaction rate for two-body recombination {A + A ^ A), 
respectively. The parameter a is the rate of reproduction or fission (A 
A + A) and r is the particle decay rate {A 0). Similarly to the pair- 
annihilation case worked out above, the additional bare terms parametrized 
by ^, b and a that we have introduced in Eq. H41|l. will be required in 
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order to carry out consistently the one-loop renormalization of the effective 
potential. Nevertheless, after this renormalization is performed, we are free 
to set the finite and renormalized values of all these couplings identically 
to zero: £^ = b = a — 0. We will so do in the final renormalized expression 
for the effective potential obtained below. Once these taken to be zero, the 
renormalized action (|41|l will correspond to the Gribov process |25) . 

Evaluating the array in |^ for this action, taking constant fields and 
then evaluating the Fourier transform yields the following matrix (note this 
is independent of both a and b) 

M{k,Lu;(t),4>)u = -2cr0 + 2A</>2 (42) 

M{k,u;;(f>,4>)i2 = -iu; + S_+Dk^ + a + r-2a4>-2X(l){l~24>){43) 

M(fc,c^; 0,0)21 = M*(A:,w; 0,0)12 (44) 

M(fc,c^; 0,0)22 = -2A(l-0)0. (45) 

Inserting this matrix into the expression for the one-loop effective potential 
Eq. ^ yields 



V(0,0) = a(l-0)00-r0(l-0)-A(l 

A f d'^k du fdctM{k,uj;<j),<j)) 



In — ' ' f46) 

2 J {2ttY J_^2ti VdetM(fc,w;0,l)/' ^ ' 

where, 

-detM(/c,w;0,0) = ^ [D/c^ _^ ^ + (T + r - 2(70 - 2A0(1 - 20)]2 

-I- 2A0(l-0)(-2cr0 + 2A02), (47) 

-AetM{k,uj;Q,l) = u"^ + [Dk^ + - a + rf . (48) 

The frequency integration in Eq. (|46|l is again performed immediately using 
Eq. (|HJ. We regulate and evaluate the remaining wavenumber integral 
directly in c? = 2 as above in Sec. (|2.1|l . The integral that is needed is 
similar in structure to that worked out in Eq. (|10|) . After some algebra, 
we find that the UV-regulated integral takes a form identical to that in Eq. 
H15|) . Put f = -\- r -\- a , then for the bare action Eq. (|41|l we have /a as 
written in Eq. lfT31) . where now the coefficients turn out to be 

ai = 2(/-2cr0-2A0(l- 20)) (49) 
02 = [/-2cr0-2A0(l-20)]2-t-2A0(l-0)(-2cr0 + 2A02) (50) 
A = 4a2 -a? = -16Act(1 - 0)00 -f 16A2(1 - 0)002. (51) 

At this point we have the following regulated one loop effective potential 
for the Gribov process (where F — a + r) 

V(0,0) = cr(l - 0)00 -r0(l-0)-A(l- 0)002 +60-^^0^ -ho 
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+ l(A + 2a?)-^V^}+0(A-2). (52) 

This expression is given in terms of bare (unrenormalized) parameters 
{a,r,X,^,b,a) and the dependence on the UV cutoff A is exphcit. As re- 
marked above, the terms proportional to bcj), ^(jxj) and the constant a will 
be needed for carrying out the renormalization program. Note that the 
regulated effective potential Eq. (|52|l is independent of the arbitrary scale 

/X. 

We now proceed to renormalize, i.e., absorb the divergences into the 
bare parameters . Specifically, we write 

a = cT{^l) + Aki + OiA"^) (53) 

r = r{p)+AK2 + 0{A^) (54) 

A ^ \{p) + Ak3 + 0{A^) (55) 

e = ^ifi) + Ak4 + 0{A^) (56) 

b = b{^i) + Ak5 + 0{A^) (57) 

a = a{fi)+AK + 0{A^). (58) 

The divergent coefficients Ki of the counterterms are chosen to cancel off 
the divergences in the one loop regulated contribution Eq. (|52|) . Inserting 
Eq. H53I58|I into Eq. (|52) yields the solutions 



47rL> ^2 



'^i = .-n i^t;2 (59) 



K2 = «4 = -Ma2 (60) 

47r 



K3 = -T^ln— 61) 



K5 = — K = — — A . (62) 
An 

After this step, we obtain the finite and renormalized one-loop effective 
potential for the Gribov process in d = 2: 

Vicf^J) = a(l-^)0(/)-r0(l-0)-A(l-<^)^02 



{ [2A.(1 - m 2A^(1 - 4>W] X In (^^9^) 
+ }_^A + 2al)-^^} +0(2 loops). (63) 
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The parameters appearing in this expression (and within ai , 02 and A) are 
understood here to be the finite scale-dependent ones, that is, a{fi),r{fj,) 
and X{fJ.). Their exphcit scale-dependence is determined below. As we 
mentioned above, we have set the finite renormalized values of ^, b and a 
identically to zero in this final expression for the potential. 



3.2 The one-loop RGEs 

Since V{(j), (j)) does not depend on the arbitrary scale /i, inserting Eq. H63() 
into Eq. H23|) will immediately yield the one-loop renormalization group 
equations in d = 2: 

= Q + 0{A^) (65) 

djjL 

As we saw above in Sec. (I2.2f) . the one-loop RGE's follow simply from 
differentiating the explicit /^-dependence in the logarithm term. The above 
three RGEs result from the vanishing of the individual coefficients of the 
independent field polynomials (1 — (f>)4>4i, {1 — (f>)<j) and (1 — 0)^0^, respec- 
tively. These equations can be solved immediately and give 

a{^l) = a„{-t^)^ exp{0{A^)) (67) 
Mo 

r{fi) = ro + OiA^) (68) 

A(m) = ^o(l-^ln^ + 0(-42))"'. (69) 

If we retain the dependence of Eq. on the other three parameters ^, b 
and a and work out the consequences implied by Eq. we would find 

that they satisfy RGEs similar to that for r in Eq. (|65|l . We conclude that 
neither ^, b nor a run with scale (at one- loop order). So, we are justified 
in setting ^(/i) = &(/i) = a{^i) = after the renormalization procedure is 
properly carried out. 

Although the RGE's in Eq. (|64I66(I have been calculated directly in 
d = 2 dimensions, as in Sec. (|2.2(l simple dimensional analysis can be 
exploited in order to infer the approximate form that these equations will 
take in other dimensions. Dimensional analysis of the action Eq. H41fl 
yields (and after absorbing the diffusion constant I? by a re-scaling of the 
time) 

[t] = [<J/D]=p\ 
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[^,|D] = [t/D] = [C/C] = [X/D] = v"-'. (70) 

From these relations, we are led to introduce the two dimensionless param- 
eters (e = 2 — d) 

We now write Eq. and Eq. so they are dimensionally correct, 

albeit only approximate in any dimension d (we herewith set ^ = 1): 

do- Kd , , 

LL— = — — a\ 72) 

where Kd is defined in Eq. H32|l . Inserting s and g from Eq. I|71|) into Eq. 
H72() and Eq. I|73|l yields the following (approximate if d 7^ 2) one-loop beta 
and flow functions 

m = ti^ = -eg + Kd9' (74) 
lis) = /-^ = -2 + Kd9. (75) 

3.3 Effective equation of motion 

Just as for the case of pair annihilation treated above, we can recover the 
effective action from knowledge of the effective potential and obtain the 
noise-corrected equations of motion for the Gribov process. Inserting Eq. 

into Eq. (|^ . we again find that cj) = 1 solves Eq. l(T7|l identically 
while the equation of motion for Eq. H38|) -the noise-averaged particle 
density- works out to be 

BV 

{dt-Dy^)<P = _ (76) 



d(j) 

-^-(Aa^-AV)ln C^-^V'^ )+0(^-) 



(Tip — r(j) — 



= F^. 

In Fig. 121 we show an example of the force experienced by the field for 
a case with r < a and equilibrium density cj) = (cr — r) /A. Setting ^ = 
on the right hand side of Eq. H76|l yields the so-called rate equation which 
one might write down as a first approximation to the equation of motion 
for the density. If we let n{t) stand for the spatial average of 0, then 
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we see that the only stationary solutions of Eq. (|76|l for zero noise are 
n = 0, corresponding to the inactive state with no particles remaining, and 
n = the active state, provided of course that a > r. Notice that for 
(j) < {a — r) /2X the effective renormalized force is ill-defined (it becomes 
complex, due to the logarithm term in the renormalized expression). But 
suppose that we arrange for the system to have initially equal rates of fission 
and decay: a = r. Then in the absence of noise, the only stationary solution 
of Eq. ifTHjl at tree- level is the inactive state n = (see Fig 2-right, tree- 
level force). However, the situation changes once the one-loop correction 
is included. Even when u = r, there can be noise-induced active states, 
corresponding to finite positive values of n, which are given implicitly by 
the transcendental equation 



(<T — An) / 2Xn \ 



See the transitions to case 1 {cr/D ~ 13.8) and to case 2 {cr/D — 15.5). 
Since the physical density must be nonnegative, for Eq. lf77|) to have pos- 
itive solutions n > requires that either (i) a > Xn > or (ii) 
a < Xn < In the first case, we are in effect canceling a term of 
order unity against a term of order cr, and there is no reason why a need 
be small: we are not employing u as a perturbation expansion parameter. 
Thus the finite density solutions satisfying criteria (i) can be valid in the 
context of perturbation theory. The solutions satisfying criteria (ii) how- 
ever fall prey to the same criticism as the solution given in Eq. (|40|l . In 
latter case, we are in effect cancelling a term of order A against a term of 
order A^. 



4 Branching and annihilating random walk 

A + A^0 and A ^ {m + 1)A 

In our third and final application of these calculational methods, we cal- 
culate the renormalized one-loop effective potential for field theories of 
branching and annihilating random walks (BARW). BARW describe the 
dynamics of a single particle undergoing three basic processes: diffusion, 
an annihilation reaction A — > and a branching process A {m+l)A, 
where m is a positive integer. A field theory for these processes was con- 
structed and its asymptotic properties analyzed in j26| . where further refer- 
ences on BARW can be found. Below we deduce the corresponding RGB's 
in d = 2 dimensions for any even m and then specialize to the case m = 2 
and solve them. The associated one-loop potential and effective force are 
plotted which show qualitatively the role of small amplitude noise on the 
system. 
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Figure 2: One-loop effective and tree-level force F{<l)). Depending on the 
values of renormalized model parameters cr(/Lt), A(yu), the noise can change 
the number of zeroes in the force from the one at the origin (tree-level force) 
to two (case 1) to three (case 2). Zeroes at finite (/) signal the presence of 
active states. 



4.1 One loop effective potential 

The field theoretic action for even-m BARW we will use in the following 
calculations is adapted from and reads as follows: 

S[^,^] = J d'^xdt i:(dt + £.- D\/^'^^/j~ X{l-4>^)i;^ (78) 

m/2 m/2 

+ ^a2i(l-^2')# + ^6z^2' 

1=1 1=0 

As before, A represents the rate of pair annihilation and the (721 stand for the 
branching rates into 2Z -I- 1 particles, for 21 = 2, 4, 6, m. As pointed out 
in j26| . the original branching process for a given fixed integer m generates 
all the lower branching reactions with m — 2,m — 4, ... offspring particles 
by virtue of the fluctuations. All these additional lower order branching 
reactions must therefore be included in the bare action Eq. (|78|l in order to 
be able to carry out a consistent renormalization procedure. In addition, 
the terms proportional to ^ and to the ^2; that we have introduced into 
l|78|l are needed to cancel off specific quadratic and logarithmic divergences 
that arise in the one-loop potential. Once we renormalize this theory, we 
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are free to set their finite values to zero, and we will do so below. Going 
through the by now familiar steps in Eq. |^ and in Eq. we find the 
action Eq. H78|l yields the matrix M{k,uj; 0, 0) of the Fourier transform of 
the second order field variations, whose individual elements are given by 



m/2 

1=1 

m/2 

m/2 



(79) 



M{k,uj-(t), (f>)i 



M{k,uj;(f>, (j))2i 
M{k,uj;(f>, ^)22 



C + Dk^ + + ^^2/(1 - [2/ + iW^) 



1=1 



M*{k,UJ-(j),(f))l2 

-2A(f - 4?). 



(80) 
(81) 
(82) 



Inserting this into the expression for the one-loop effective potential Eq. 
© yields 



%/2 



-1/2 



V(0,0) = 

+ 

where 

- detM(fc,a;;( 

+2A(l-02) 
and 



-A0'(l-0') + ^a2i(l 



1=1 



1=0 



A 
2 



d'^k duj 



In 



dct M{k, uj; cj), (j)) 



{2TTf2n VdetM(fc,w; 0,1)7' 



V2 



(83) 



DP + ^ + 4A# + (72,(1 - [21 + 1]^2') 



1=1 



m/2 m/2 

2\cf>^ - Y <^2i{2l + l)(2/)02/-i^ ^ J2 6/2/(2/ - 1) 

1=1 1=0 



T21-2 



m/2 

det M(fc, cj; 0, 1) = + [Dk^ + $ - ^ (2/)(72,] ' 



(84) 



The frequency integration is again immediate, and the remaining (UV- 
regulated) integral over wavenumber fc in d = 2 has a mathematical struc- 
ture identical to that in Eq. (|10|l . although we need to simply replace the 
(u -l- 6^ term by [u -\~ F)"^ under the square-root in the second integrand. 
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where F = ^ — Y^^i 2^o'2/- The field dependent coefficients appearing in 
/a in Eq. (|10|l corresponding to even m BARW are thus given by 



m/2 



ai = + 4A# + ^ (72/(1 - [2; + l]02/)j 



(85) 



m/2 



a2 = ^ + 2A(l-<^2) 2A02„^^2;(2/ + l)(2O02'-\ 



m/2 

+ ^6;2;(2?-l) 



72/-2 



/=0 



(86) 



m/2 



A = 4a2-a? = 8A(l-</.2) 2A</.2-^CT2i(2? + l)(2O0^'"\ 



(=1 



m/2 

Y,^2i2l{2l-1) 

1=0 



T21-2 



(87) 



Examining the ultraviolet divergences in this regulated wavenumbcr inte- 
gral for arbitrary even m, reveals the same general pattern of quadratic and 
logarithmic divergences in addition to finite terms as spelled out explicitly 
in Eq. p5|l . From Eq. (|85I87|) . we immediately see that the polynomial 
field structure of these two divergences contains operators of the same order 
as those appearing in the tree-level potential Vq for arbitrary even m. The 
coefficients of these quadratic and logarithmic divergences can be written 
as 



A 



i/2 



= 4A^(/) + ^a2i[2; + l](l 



(88) 



1=1 



m/2 



1=1 

T2i-2 



m/2 

^6/2/(2^-1) 
2A2(1-0^)(/)2 

m/2 

A ^[-(20(2^ + l)<J2i + {21 + 2){2l + 3)(72i+2](l 



12/ 



i=l 



m/2 



2A6 + A ^[(2Z + 2){2l + 1)^21+2 - 2^(2/ - 1)6;]^'', (89) 



1=1 
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respectively, where we define ct„i+2 = and = 0. For all even m we 

conclude that this field theory is (at least) one-loop renormalizable. Using 
Eq. H15|l once again to separate the divergences from the finite parts in the 
logarithm allows us to write the regulated one-loop potential as follows, 

m/2 m/2 

v(0,0) - -X{i~^')cb^ + Y,<^2i{i~4>^')4><l> + Cci,^ + Y.^^i^^' 

+ -L(A + 2a?)-^V^| + 0(A-2). (90) 

We renormalize this potential by absorbing the divergences into the bare 
parameters and introduce a set of countertcrms, where I = (1,2, ...,m/2) 

X = X{n) + Akx + 0{A^) (91) 

+ Ak^,, + 0{A^) (92) 

^ = Ci^i) + Ak^ + 0{A^) (93) 

6; = i2i{p) + Ak^,, + 0{A^) (94) 

Co = M+A^^io+OiA''). (95) 

The following choice of countertcrms 
47rL> fi^ 



= _A^ln_ (96) 



X{fi)[2l{2l + l)aM - {21 + 2){2l + 3)a2,+2(M)] , A 



2 



8nD iJL 



In^ (97) 



A 2 



H = -^A^ (98) 



m/2 



(2[+lWM .2 

- ^[(2/ + 2)(2; + 1)6/+2(m) - 2/(2/ - 1)6Km)] In ^,(100) 

(recall (7^+2 ~ Cm,+2 = 0) where / = 1,2, ...,m/2, yields the renormalized 
one-loop effective potential for all even m BARW: 

m/2 

V((/.,0) = -A(/i)(l- 02)02 
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m/2 

+ < In -^^ — +— (A + 2a?) -y^ > 

8ttD\ 8 V 4L>^2 y ^ ^ 4^ J 

+ 0(2 Zoops). (101) 

Here, ai, 02 and A are given in Eq. (|85l8t)|) and Eq. (|57|l . respectively and 
the parameters A(^), cr2i{fJ.) and £,2i{p) appearing in this final expression 
are scale dependent. We determine their scale dependence below. 



4.2 The one-loop RGEs 

The renormalized potential Eq. H101() does not depend on the arbitrary 
scale As a result, inserting Eq. (|101|l into Eq. H23|l will immediately 
yield the one-loop renormalization group equations in d = 2 dimensions: 

dX 1 ,9 

^T, = 2^' ^'''^ 

= -^[2^(2/ + 1)^2/ -(2? + 2)(2/ + 3)(72,+2] (103) 
fi^ = + O{A^) (104) 

= -2^ 

= -^[2^(2/ -1)6, -(2^ + 2) (2; + 1)6^+2]. (106) 

As we saw above in Sec. H2.2II and in Sec. (|3.2|l . at one- loop order, the 
RGE's follow immediately from differentiating the explicit /i-dependence 
in the logarithm term in Eq. H101|l . The above set of RGEs results from 
picking off the coefficients of the field polynomials (1 — <P)<p'^, (1 — (p'')(j)<j>, 
for I = 1, 2, r7i/2, (fxp and of cfP'' , ior I = 0, 1, 2, to/2, respectively. 
Note that the RGE Eq. pU3|l at the highest value Z = to/2 is given by 

dam m{m + l) 

Up to this stage, we have carried through the complete calculation of the 
effective potential valid for arbitrary even to. However, it is known from 
Ref. that the most relevant of the branching reactions is actually the 
one with smallest to, namely to = 2. The branching process with m — 2 
therefore will describe the entire universality class of BARW with even 
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offspring. For this reason, we now specialize the remainder of this section 
to this case. 

The explicit solutions of the m — 2 RGB's Eq. (|102llUt)|l are as follows: 

3A>.Q 

<J2it^) = a2{0)(^)"'" exp{OiA^)) (109) 

c(/i) = m + o{A^) (110) 

?o(M) = U0)-^^lnJ^ + O{A^) (111) 

e2(M) = e2(0)(^)'^exp(O(^2))^ (^^2) 

According to these solutions, we can now set the renormalized values of 
^(a*);Co(/^) and £,2{fJ-) all identically equal to zero at all scales by simply 
choosing their values at the arbitrary initial scale fio to be zero: ^(0) = 
^o(O) = ^2(0) = 0. Our final renormalized potential depends only on the 
rate of pair annihilation and the branching rate. The potential and effective 
force will be calculated with this choice (see Sec (|4.3f) below). 

The RGB's Bq. (|102I106|I have been calculated directly in d = 2, which 
coincides the upper critical dimension for BARW. As before, we can deduce 
crude and approximate forms that these RGB's must take in any dimension 
by combining simple dimensional analysis with these d — 2 RGB's. We take 
Bq. HlU2fl and Bq. (|107|) for m = 2 and write them as follows: 

^ (113) 

a/i 2'kD^'^ 

da2 3 ^ (114) 

dfi 2'nDpL'^ 

Then introduce the dimensionless parameters 

into Bq. H113I114|I . This yields the approximate one-loop beta and zeta 
functions for m — 2 BARW in dimension d: 
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For d = 2, these reduce to 




2 + 3g, 



(119) 



(120) 



and reproduce the d — 2 one-loop beta and zeta functions obtained in PS| 
by other means. 

4.3 Effective equation of motion 

We can calculate the effective force and equation of motion for the physical 
density 4> starting from Eq. H3t)|l (note there is no wavefunction renormal- 
ization for BARW ^^)- Once again, we find that = 1 solves identically 
its equation of motion Eq. ()37fl , while evolution of is determined by by 
Eq. which when written out yields 



(j) with zeroes at cf) — aijX. The one-loop force is ill defined when the 
argument of the logarithm goes negative. Here, if we set directly (72 = in 
the force in Eq. H121|) . the one- loop correction does not induce any active 
states that can be considered valid within perturbation theory, for the same 
reasons as spelled out in the case of pair-annihilation. 



5 Discussion 

In this paper we have taken the well established effective action and effec- 
tive potential formalism from the quantum field theory domain and have 
adapted and applied it to classical field theories of the Doi and Peliti type of 
reacting and diffusing particles. These are field-theoretic representations of 
the corresponding diffusion-limited reaction problems derived from classi- 
cal master equations. Apart from taking the continuum limit, the resultant 
actions require no further approximation and the noise is included exactly. 
The main result of our application is the formula for the one-loop effective 
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Li. 




1) 

Figure 3: One-loop effective and tree-level force F{(j)). The zero at the ori- 
gin corresponds to the inactive state, while the zero at finite (j) corresponds 
to an active state. 



potential which holds for any multiple A^-species reaction-diffusion system 
subject to internal particle density fluctuations. The main goal of this pa- 
per is to illustrate the relative ease with which this effective potential can 
be analytically calculated in closed form and consistently renormalized for 
a wide variety of interesting single {N = 1) species models, of increasing 
complexity, from the simplest pair annihilation case to the more intricate 
case of branching and annihilating random walks with arbitrary even num- 
ber m of offspring. We believe the most interesting specific application of 
the effective potential formalism in these problems is the ease with which 
one can deduce the associated one-loop renormalization group equations 
(RGE) and the effective equations of motion from it. 

We have calculated the effective potential in fixed space dimension 
d — 2, corresponding to the upper critical dimension for all the models 
surveyed here. Once the frequency integration is carried out in Eq. we 
have seen that the remaining integral over the wavevector always reduces 
to the same small set of simple elementary functions (polynomials and log- 
arithms), irrespective of the microscopic details of the specific model. This 
feature is what allows one to obtain potentials and the RGE's for a wide 
variety of different models with little additional effort. In so far as this is 
a consequence of the one-loop approximation, this does represent a major 
technical advantage over the standard Feynman diagram approach, since 
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in the latter approach, each individual model requires the handling of a 
distinct set of diagrams with its attendant combinatorics to work out. Of 
course, when one loop order is not sufficient, to go to higher loops, one must 
revert back to diagrams to compute the higher order corrections to the po- 
tential, and this can be done While individual Feynman diagrams can 
be analytically continued to arbitrary dimension, this appears to be more 
difficult for the potential, representing, as it does, an infinite summation 
of diagrams |221- Nevertheless, the analytic continuation of the one- loop 
potential to arbitrary dimensions may be possible, and this point clearly 
deserves further work. This would open up the very interesting possibility 
to derive the associated RGB's for a wide variety of reaction-diffusion mod- 
els in a straightforward manner for any dimension d in a universal way, at 
one-loop level, without the need to handle individual Feynman diagrams. 
It would be interesting to test out these analytic methods on two-species 
{N = 2) reaction-diffusion systems such as the A + B ^ reaction [77] . 
The expression for an = 2 species potential requires working out the 
determinant of 4 x 4 arrays, followed by integrations over frequency and 
wavenumber. Our conviction is that further developments of these poten- 
tial methods can promote significant advances in noncquilibrium statistical 
mecanics and in physical chemistry. The simple examples treated here serve 
as a testing ground for this technique, and we believe we have amply demon- 
strated the usefulness of the approach and the relative ease in carrying out 
the specific steps leading to the renormalization group equations. 



A Effective potential 

In this appendix we derive the one-loop effective potential Eq. (|3J) for 
Doi-Peliti type stochastic field theories. Following standard procedures 
ISIEIj one starts from the microscopic master equation corresponding to the 
specific chemical reactions to be studied and then passes to a field-theoretic 
solution of this master equation. In the final step, one obtains a path 
integral representation of the time-evolution operator (or "hamiltonian" ) 
of the general form 

J nf (PV.PV^,) exp ( - S[A, ^^]) , (122) 

for A^-species of particles ipi and N conjugate or "response" fields -ipi for 
i = (1, 2, ... , N), and where S denotes the unrenormalized (bare) action. 
Introduce the loop counting parameter ^ = 1 in the exponential in Eq. 
(|122|l as follows 

S[i;^,i;^]^ S[iP,,^ij^]/A. (123) 
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Fluctuations are taken into account via the loop expansion in powers of A. 
Since the loop expansion corresponds to an expansion in a parameter that 
multiplies the entire action, it is unaffected by shifts of fields. 

We next adapt the basic steps that take us from an action S to the 
effective action and effective potential P □ [H UHl OOl EH EH ■ First, we 
introduce a set of arbitrary source functions Ji and Jj for both the ipi and 
tpi by adding the following term to the action 

^ S[^„^,] + jd^'x J dt Y,{X,^, + J.^.), (124) 

this linear shift leads to a source dependence in the path integral Eq. (|122|) 
which we denote by Z[J, J]. Next, introduce the functional for connected 
correlation functions W[J, J] and its Legendre transform, the effective ac- 
tion T[(f), (f)] 

W[J,J] = A\nZ[jJ], (125) 
r[0,0] = -W[J.,J]+ I d'^x [ dtyU,cp, + J,^,), (126) 



i=l 



where 



<Pr[J,J\ = — Yl , 0,;[J, JJ = Jj . (127) 

OJi oJi 

The fields 0^ and Ji (and and Ji) make up the Legendre transform pairs. 
From Eq. (|122I124I125() . and Eq. lfn7|l . we see that the fields are 
the fluctuation averages of the original degrees of freedom, but evaluated 
in the presence of the source terms, that is, 

0,[J, J] = (^,(x, t))jj 0,[J, J] = (Mx, t))jj , (128) 

where the angular brackets denote the average over the fluctuations. These 
fluctuation-averaged fields are the solutions of the effective equations of 
motion 



_ _ _ = J„ (129) 
0[J„/]>[J,J] 



in the presence of the source terms, a fact which follows immediately from 
differentiating Eq. (|126|l as indicated. Most importantly, once we set the 
external sources to zero in Eq. H129(l . J, J = 0, we then obtain the effective 
equations of motion for the stochastic-averaged fields. The solutions are 
the stationary points of the effective action: 

(3M) =„, =0. (130) 
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We make use of Eq. (|130|) in the paper in specific reaction diffusion models 
and explicitly obtain the fluctuation-modified effective equations of motion 
satisfied by the noise-averaged particle densities at one-loop order. 

To complete the derivation of the one-loop effective potential (and sup- 
pressing writing out the full dependence on indices, and spatial and tem- 
poral arguments) we henceforth let {(j), (f>) denote the set of the source-free 
fluctuation averages of the 2N scalar fields {(i/'i)(o,o) i (V'i)(o,o)}fc=i- Then 
the effective action for the field theory described by Eq. (|122|l is given by 

1201 

m 4>] = S[cj,, + AT, + OiA^), (131) 
where the one-loop contribution is explicitly calculated to be 

Ti ^] = iTr[Ln S^^^ (</., ^; x, y) - Ln S^^'> (0, 1; x, y)], (132) 

and Tr and Ln are functional traces and logarithms (these operate on both 
discrete and continuous indices). The array S'(^) of second order functional 
derivatives of the action with respect to the fields is expressed in 2 x 2 
block array form in Eq. where i,j = (1,2,..., N). The action S[(j), (f>] 
is simply the action Eq. H123() after replacing ipi and tpi by 4>i and 0i, 
respectively. The effective potential is the effective action Eq. H131|) evalu- 
ated on constant field configurations, and after dividing out by the overall 
volume of space-time VT: 

IV, V\ yrp VT VT ^ ' 

= V^[<P,$\+AV^[ct),^]+0{A^). (133) 

For constant fields 4'i4' = const., we pass to Fourier space and define a 
2N X 2N matrix M via the Fourier transform (FT) of Eq. (Q) as defined 
in Eq. (0). Since this expression is diagonal in momentum and frequency, 
we can straightforwardly carry out the sequence of operations indicated in 
Eq. ifT^ lis,: 

TrLn[M(fc, ui; 0, 0) x (27r)'^+M''(fc - p)5{oj - 9)] 
= Tr(27r)'^+i(5'*(A; - p)6{uj - n) In M{k, oj, cf), 0) 

— / —{27rf+'S^ik^p)S{oj-n)\k=p.^^ntrlnM{k,co,^,4>) 

[ZTT ) J Zn 

d'x Jdt J ^ I ^lndetM(fc,c^,0,0), (134) 

where In and tr are the ordinary logarithm and trace, and we have used 
the identity trlnM = IndetM. Evaluating Eq. (|131|l on constant config- 
urations and dividing through by the volume of space-time and using Eq. 
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H1I2I132I133|) and Eq. p34l) , yields the compact expression for the one- loop 
effective potential Eq. Q. 
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